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By Peter Brockwell, 1 Erdenebaatar Chadraa 
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A family of continuous-time generalized autoregressive condition- 
ally heteroscedastic processes, generalizing the COGARCH(l, 1) pro- 
cess of Kliippelberg, Lindner and Mailer [J. Appl. Probab. 41 (2004) 
601-622], is introduced and studied. The resulting COGARCH(p,q) 
processes, q > p > 1, exhibit many of the characteristic features of 
observed financial time series, while their corresponding volatility 
and squared increment processes display a broader range of autocor- 
relation structures than those of the COGARCH(l,l) process. We 
establish sufficient conditions for the existence of a strictly station- 
ary nonnegative solution of the equations for the volatility process 
and, under conditions which ensure the finiteness of the required mo- 
ments, determine the autocorrelation functions of both the volatility 
and the squared increment processes. The volatility process is found 
to have the autocorrelation function of a continuous-time autoregres- 
sive moving average process. 

1. Introduction. In financial econometrics, discrete-time GARCH (gen- 
eralized autoregressive conditionally heteroscedastic) processes are widely 
used to model the returns at regular intervals on stocks, currency invest- 
ments and other assets. Specifically, a GARCH process (£ n )neN typically 
represents the increments lnP n — lnP n _i of the logarithms of the asset price 
at times 1,2,3,... . These models capture many of the so-called stylized 
features of such data, for example, tail heaviness, volatility clustering and 
dependence without correlation. For GARCH processes with finite fourth 
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moments, the autocorrelation functions of both the squared process and 
the associated volatility process are those of autoregressive moving aver- 
age (ARMA) processes. The squared GARCH(1,1) process, for example, 
has the autocorrelation function of an ARMA(1,1) process and the corre- 
sponding volatility has the autocorrelation function of an AR(1) process. 
Higher-order GARCH(p, q) processes were introduced to allow for the possi- 
bility of a broader range of autocorrelations for the volatility and the squared 
increment processes. 

Various attempts have been made to capture the stylized features of finan- 
cial time series using continuous-time models. The interest in continuous- 
time models stems from their use in modeling irregularly spaced data, their 
use in financial applications such as option pricing and the current wide- 
spread availability of high-frequency data. In continuous time it is natural 
to model the logarithm of the asset price itself, that is, Gt = lnPt, rather 
than its increments as in discrete time. 

Notable among these attempts is the GARCH diffusion approximation of 
Nelson [23]. (See also [12] and [13].) Although the GARCH process is driven 
by a single noise sequence, the diffusion limit is driven by two independent 
Brownian motions (W t )t>o and (W^)t>o- For example, the GARCH(1, 1) 
diffusion limit satisfies 

dG t = a t dW t {1 \ 

(1.1) 

do% = 0(i-o$) + po%dwl 2) , t>0. 

The behavior of this diffusion limit is therefore rather different from that 
of the GARCH process itself since the volatility process (a 2 )t>o evolves 

independently of the driving process {W^)t>Q in the first of the equations 
(1.1). 

Another approach is via the stochastic volatility model of Barndorff- 
Nielsen and Shephard [3, 4] in which the volatility process a 2 is an Ornstein- 
Uhlenbeck (O-U) process driven by a nondecreasing Levy process and G 
satisfies an equation of the form dGt = pdt + at dWt , where If is a Brow- 
nian motion independent of the Levy process. The autocorrelation function 
of the Levy-driven O-U volatility process has the form p(h) = exp(— c\h\) 
for some c > 0, but this class can be extended by specifying the volatility 
to be a superposition of O-U processes as in [2] or a Levy-driven CARMA 
(continuous-time ARMA) process as in [10]. As in Nelson's diffusion limit, 
the process G is again driven by two independent noise processes and the 
volatility process a 2 evolves independently of the process W in the equation 
for G. 

A continuous-time analog of the GARCH (1,1) process, denoted 
COGARCH(l, 1), was recently constructed and studied by Kliippelberg, 
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Lindner and Mailer [19] . Their construction is based on taking a limit of an 
explicit representation of the discrete-time GARCH(1, 1) process to obtain a 
continuous-time analog. Since no such representation exists for higher-order 
discrete-time GARCH processes, a different approach is needed to construct 
higher-order continuous-time analogs. In this paper we do this by specifying 
a system of Levy-driven stochastic differential equations for the processes 
G and a 2 . If the volatility process a 1 is strictly stationary we refer to the 
process G as a COGARCH(p, q) process. In the special case p = q = 1 we re- 
cover the COGARCH(l, 1) process of Kluppelberg, Lindner and Mailer [19]. 
In general we obtain a class of processes G with uncorrelated increments 
but for which the corresponding volatility and squared increment processes 
exhibit a broad range of autocorrelation functions. The volatility process 
has the autocorrelation function of a continuous-time ARMA process. 

The construction of the COGARCH(l,l) process due to Kluppelberg, 
Lindner and Mailer [19] starts from the defining equations of the discrete- 
time GARCH(1,1) process (£ n )neN > 



(1.2) 



0~n 



«0 + aiCn-i + Picrl-i, 



n e N , 



where ao,«i and (3± are all strictly positive, and (e n )neN is a sequence of 
i.i.d. (independent and identically distributed) random variables with mean 
and variance 1. The recursions (1.2) can be solved to give 
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cr + a / exp 



W 

+a 1 e 2 j ) 



j 



ds 



x exp 



n-l 

^log(/?i +aiSj) 



where [s\ denotes the integer part of s £ M. The C0GARCH(1, 1) equations 
are then obtained by replacing the driving noise sequence (e n )neN by the 
jumps (AL 4 = L t — L t ~)t>o of a Levy process. More precisely, observing 
that Ep^log^i + aie]) = n\ogf3 l + EpoM 1 + (ai/A)^ 2 ) for /3i > 
and writing r] for — log/5i, ujq for qq and lo\ for a\, leads to the equations 



(1.3) 

(1.4) 

where 
(1.5) 



dG t = a t dL t 



o\ = ^<7q + ojq J e Xs ds^j e 



t > 0, G = 0, 
t >0, 



X t :=rjt- log(l + cuie"(AL s ) 2 ). 

0<s<t 
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Here, loq > 0,uj\ > 0,7/ > and Og is independent of (L t ) t > . The 
C0GARCH(1,1) process is the solution G of these equations and, under 
specified conditions on the coefficients and the distribution of Og, the volatil- 
ity process a 2 is strictly stationary and G has stationary increments. 

The COGARCH(l,l) process with stationary volatility has been shown 
to have many of the features of the discrete time GARCH(1, 1) process. As 
shown in [19, 20], the COGARCH(l, 1) process has uncorrelated increments, 
while the autocorrelation functions of the volatility a 2 and of the squared in- 
crements of G decay exponentially. Furthermore, the COGARCH(l, 1) pro- 
cess has heavy tails and volatility clusters at high levels; see [14] and [20]. 
For an overview of the extremes of stochastic volatility models, see [14] and 
[21]. Also, observe that many of the features of the COGARCH(l, 1) process 
can be obtained in a more general setting, as in [21]. 

In the next section we specify a system of stochastic differential equations 
for the COGARCH(p, q) process G and its associated volatility process, 
which we shall denote by V. This is directly motivated by the corresponding 
structure of the discrete-time GARCH(p, q) process. We then show that 
the solution of these equations coincides with that of the COGARCH(l,l) 
equations if p = q = 1. Notation and definitions used throughout the paper 
are given at the end of Section 2. 

In Section 3 we give sufficient conditions for the existence of a strictly 
stationary volatility process. In the COGARCH(l, 1) case, these are exactly 
the necessary and sufficient conditions obtained by Kluppelberg, Lindner 
and Mailer [19, 20]. More detailed results are given in the special case when 
the driving Levy process is compound Poisson. The proofs rely on the fact 
that the state vector of the COGARCH(p, q) process, sampled at uniformly 
spaced discrete times, satisfies a multivariate random recurrence equation. 

In Section 4 we focus on the autocorrelation structure of the stationary 
volatility process. Just as the discrete-time GARCH volatility process has 
the autocorrelation function of an ARMA process, the COGARCH volatility 
process has the autocorrelation function of a CARMA process. 

Section 5 deals with conditions which ensure positivity of the volatility, 
while the autocorrelation structure of the squared increments of the COGA- 
RCH process itself is obtained in Section 6. The results are illustrated with 
simulations in Section 7. So as not to disturb the flow of the arguments, 
proofs of the results are postponed to Sections 8-11. 

2. The COGARCH(p, q) equations. Let (e n )neN De an i-i-d- sequence 
of random variables. Let p, q > 0. Then the GARCH(j>, q) process (Cn)neN 
is defined by the equations 

(2.1) a 2 =ao + ai£n-i H h a p f£_ p 
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+ Pio%-! H \-P q al_ q , n>s, 

where s := max(p, q), tTg, . . . , cr^_ 1 are i.i.d. and independent of the i.i.d. 
sequence (e n ) n > s , and £ n = G n+ \ — G n represents the increment at time 
n of the log asset price process (Gn)neN - Note that the continuous-time 
GARCH process will be a model for (Gt)t>o an d not for its increments as 
in discrete time. 

Equation (2.1) shows that the volatility process (c>"n)neNo can De viewed 
as a "self-exciting" ARMA(q,p — 1) process driven by the noise sequence 
(o"n-i e n-i)neN- Motivated by this observation, we will define a continuous- 
time GARCH model for the log asset price process (Gt)t>o of order (p, q) 
by 

dG t = a t dL u t>0, G = 0, 

where (&t)t>0 is a CARMA(<7,p — 1) process driven by a suitable replacement 
for the discrete time driving noise sequence (&n-i £ n-i)neN- 

The state-space representation of a Levy-driven CARMA((/,p— 1) process 
(^t)t>o with driving Levy process L, location parameter c, moving average 
coefficients a.\, . . . , a p , autoregressive coefficients f3i,...,/3 q and q > p is (see 
[8]) 

ipt = c + a'Ct, 
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where a' = [a± , . . . , a q ] , aj := for j > p and the coefficient matrix in the last 
equation is —/3i if 9 = 1. (The CARMA(g,p — 1) process (V ; t)t>o is a strictly 
stationary solution of these equations, which exists under conditions found 
in [9].) To obtain a continuous-time analog of the equation (2.1), we suppose 
that the volatility process {aJ)t>o has the state-space representation of a 
CARMA(<7,p — 1) process in which the driving Levy process (Lt) is replaced 
by a continuous-time analog of the driving process (Cn-i £ n-i)neN in (2.1). 

The increments of the driving process in continuous time should corre- 
spond to the increments of the discrete-time process: 

n—1 n—1 

i=0 i=o 

We therefore replace the innovations e n by the jumps ALt of a Levy process 
(Lt)t>o to obtain the continuous-time analog 

R t := Y, <^-(AL,) 2 , *>0. 

0<s<t 
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If L has no Gaussian part [i.e., r| = in (2.2) below], we recognize R as the 
quadratic covariation of G, that is, 

R t= E ^ 2 -(AL s ) 2 = [ t o 2 a _d[L,L] a = [G,G\ t . 

If L has a Gaussian part, then still J^OKsKti^-^s) 2 = [L,L]^, the discrete 
part of the quadratic covariation, and we have in general 

Rt = / of- , that is, dR t = a 2 t _ d[L, L][ d) . 

Jo 

Recall that for a Levy process L = (Lj)t>o the characteristic function 
E(e tSLt ),6 G R, can be written in the form 

, jyC i6L t , (J^.a -2? 

(2.2) 



E{e w ^) = eMt[i lL e-ri 



+ / (e lOx -l-i0xl lxl<1 )dv L (x 

JR 

The constants 7^ S R, t\ > and the measure ul on R form the charac- 
teristic triplet of L. As usual, the Levy measure is required to satisfy 
J K min(l,a; 2 ) di>i,{x) < 00. For more information on Levy processes, we refer 
to the books by Applebaum [1], Bertoin [5] or Sato [25]. 

The COGARCH(p, q) equations will now be obtained by specifying that 
the volatility process V (= a 2 ) should satisfy continuous-time ARMA equa- 
tions driven by the process R defined above. Provided V is nonnegative 
almost surely (conditions for which are given in Section 5), we can define a 
process G by the equations Go = and dGt = s/VtdL t . Under conditions en- 
suring that V is also strictly stationary, we refer to G as a COGARCH(p, q) 
process. As we shall see, when p = q = 1, the solution of the COGARCH 
equations coincides with that of the COGARCH(l, 1) equations (1.3)-(1.5) 
of [19] . [The parameters j3\ , . . . , f3 q and a,\ , . . . , a p in the following definition 
should not be confused with the parameters denoted by the same symbols 
in the defining equation (2.1) of the discrete-time GARCH process.] 



Definition 2.1 [The COGARCH(p,g 
gers such that q > p > 1, ao > 0, a±, . . . , a % 



equations] . If p and q are inte- 

G R, /3,€R, op ^0,0,^0 



and afp-t-i 
and e by 



OLr, 



B 



0, we define the (q x q) matrix B and the vectors a 








-Pa 



1 




-/V 




1 



-/V 






1 

-A. 



CONTINUOUS-TIME GARCH PROCESSES 



7 







"0" 











e = 




a q -i 











1 



with B := —Pi if q = 1. Then if L = (Lt)t>o is a Levy process with nontrivial 
Levy measure, we define the (left-continuous) volatility process V = (Vt)t>o 
with parameters B, a, ao and driving Levy process L by 

V t = a + a'Yt_ , t > 0, Vb = a + a'Y , 

where the state process Y = (Y t )t>o is the unique cadlag solution of the 
stochastic differential equation 

(2.3) dY t = BY t _dt + e(a + a'Y t ^)d[L,L] { f\ t>0, 

with initial value Yo, independent of the driving Levy process (Lt)t>o- If 
the process (Vt)t>o is strictly stationary and nonnegative almost surely, we 
say that G = (G t )t>o, given by 

dG t = Vv t dL t , t>0, G = 0, 

is a COGARCH(p, q) process with parameters B, a, ao and driving Levy 
process L. 

That there is in fact a unique solution of (2.3) for any starting random vec- 
tor Yo follows from standard theorems on stochastic differential equations 
(e.g., [24], Chapter V, Theorem 7). The stochastic integrals are interpreted 
with respect to the filtration F = (J-t)t>o, which is defined to be the smallest 
right-continuous filtration such that J-q contains all the P-null sets of 
(L t )t>o is adapted and Yo is .Fo-nieasurable. 

Without restrictions on ao,a and B, the process V is not necessarily non- 
negative, in which case G is undefined. Conditions which ensure that V is 
nonnegative will be discussed in Section 5. In particular, it will be shown 
that if &'e Bt e > for all t > and Yo is such that V is strictly stationary, 
then V is nonnegative with probability 1. Even if V takes negative values, 
however, the process is of some interest in its own right and many of our 
results for V are valid without the nonnegativity restriction. 

Conditions for stationarity of V are discussed in Section 3. 

We next show that if p = q = 1, the solution of the COGARCH equa- 
tions in Definition 2.1 coincides with the solution of the COGARCH(l,l) 
equations of Kliippelberg, Lindner and Mailer [19]. 

Theorem 2.2. Suppose that p = q = I, and that ao,ai and P are all 
strictly positive. Then the processes (Gt)t>o and (Vt)t>o of Definition 2.1 
are, respectively, the processes (Gt)t>o and (of )t>o defined by (1.3)-(1.5), 
with parameters ojq = a^Pi, uji = ct\e~^ and rj = P\. 
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Proof. From 

dY t = -foYt dt + V t d[L, L][ d) and V t+ = a + a x Y t 
it follows that 

dVt+ = ot\ dYt = —aiPi— — — dt + aiVtd[L,L]{ d \ 

and hence that 

V t+ = a Pit-Pi fv s ds + ai V V S (AL S ) 2 +V . 
Jo nrt. 



However, this equation is also satisfied by the volatility process (of )t>o 
of (1.4) when ujq = ao/?i, i\ = Pi and u>\ = a\e~^ 1 , as shown in Proposi- 
tion 3.2 of [19], and uniqueness of the solution gives the claim. □ 

We conclude this section with a few definitions and some notation which 
will be used throughout the paper. 

Definition 2.3. Let a and B be as in Definition 2.1. Then the charac- 
teristic polynomials associated with a and B are given by 

a(z) :=ai + a.2Z + • • • + a p z v i , z & C, 
b{z) := + fcz"- 1 + ■ ■ ■ + f3 qi zeC. 

The eigenvalues of the matrix B (which are exactly the zeroes of b) will be 
denoted by Ai, . . . , X q and assumed to be ordered in such a way that 

KA 9 < KA g _i < • • • < KAi 

(where 3?Aj denotes the real part of Aj). Furthermore, define 

A := X(B) :=5RAi. 

For the rest of the paper, convergence in probability will be denoted by 
, uniform convergence on compacts in probability by and equality 
in distribution by "=". For i£Kwe shall write log + (x) for log(max{l, x}). 
The transpose of a column vector c£C will be denoted by c' . If || • || is 
a vector norm in C 9 , then the natural matrix norm of the (q x g) matrix 

C is defined as ||C|| = sup ceC9 \j } ^prjp Correspondingly, for r G [l,oo], we 
denote by || • || r both the vector L r -norm and the associated natural matrix 
norm. Recall that the natural matrix norms of the L 1 ,^ and L°° vector 
norms are the column-sum norm, the spectral norm and the row-sum norm, 
respectively. 
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The (q x q) identity matrix will be denoted by I q or simply I, and the 
canonical vector (0, . . . , 0, 1, 0, . . . , 0)', with ith component equal to 1, will be 
denoted by e^. For e q we simply write e. By diag(Ai, . . . , X q ) we mean the 
diagonal (q x q) matrix with these entries on the diagonal. The Kronecker 
product of two {q x q) matrices A and B will be denoted by A (g) B, and 
by vec(A) we denote the column vector in C q which arises from A by 
stacking the columns of A in a vector (starting with the first column). For 
the properties of the Kronecker product, we refer to Liitkepohl [22]. 

3. Stationarity conditions. In this section we consider conditions under 
which the volatility process (Vt)t>o specified in Definition 2.1 is strictly 
stationary. The parameters B, a and ao, and the state process (Y()t>o are 
as specified in Definition 2.1. The condition (3.2) established in Theorem 3.1 
below is necessary and sufficient for stationarity in the special case p = q = 1. 
For larger values of p and q it is sufficient only, but not unduly restrictive, 
since there is a rich class of models satisfying the condition. Without serious 
loss of generality we shall assume that the matrix B can be diagonalized. 
Since the only eigenvectors corresponding to the eigenvalue Aj are constant 
multiples of [1, Aj, A?, . . . , A?" 1 ]', this is equivalent to the assumption that 
the eigenvalues of B are distinct. Let S be a matrix such that S~ 1 BS is a 
diagonal matrix, for example, 

- 1 ... i - 

Al • • • \ q 

(3.1) S = . 

-A? 11 ••• A* 11 . 
[For this particular choice, S'^^BS = diag(Ai, . . . , X q ).] 

Theorem 3.1. Let (Y t ) t >o be the state process of the COGARCH(p,g) 
process with parameters B, a and ao- Suppose that all the eigenvalues of B 
are distinct. Let L be a Levy process with nontrivial Levy measure ul and 
suppose there is some r£ [l,oo] such that 

(3.2) f log(l + HS-WSH^ 2 ) du L (y) < -A = -X(B) 

for some matrix S such that S~ 1 BS is diagonal. Then Yt converges in 
distribution to a finite random variable Yqo, as t — > oo. It follows that if 

Yo = Yqo, then (Yt)t>o and (Vt)t>o are strictly stationary. 

Remark 3.2. (a) If (y t ) t >o is the volatility of a COGARCH(l, 1) process 
with parameters B = —j3\ < 0, ao > and a\ > 0, then ||S ,_1 ea'S'|| r = u\ 
and, as already indicated, the condition (3.2) is necessary and sufficient for 
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the existence of a strictly stationary COGARCH(l, 1) volatility process. (See 
[19], Theorem 3.1.) 

(b) For general q > 2, the quantity ||S'~ 1 ea5'|| r depends on the specific 
choice of S and on r. Observe that it is sufficient to find some S and some 
r such that (3.2) holds. 

The proof of Theorem 3.1 will make heavy use of the general theory of 
multivariate random recurrence equations, as discussed by Bougerol and 
Picard [7], Kesten [18] and Brandt [6] (in the one-dimensional case). The 
COGARCH state vector satisfies such a multivariate random recurrence 
equation, as indicated in the following theorem. 

Theorem 3.3. Let (Y t ) t >o be the state process of the COGARCH(p, q) 
process with parameters B, a and ao, and driving Levy process L. Then 
there exists a family (Js,ti^-s,t)o<s<t °f random {q x q) matrices J s j and 
random vectors K. s t i n ^ 9 such that 

(3.3) Y t = J Sjt Y s + K M , 0<s<t. 

Furthermore, the distribution of (J s t,K Si t) depends only ont — s, (J slj t 1 ,K Sl ^) 
and ( J S2; t 2 , K S2jt2 ) are independent for < s\ < t\ < S2 < t^; and for < s < 
u<t, 

(3.4) J S)t = J u ,tJs,u- 

If additionally the conditions of Theorem 3.1 hold, then the distribution of 
the vector Yoo is, for any h > 0, the unique solution of the random fixed 
point equation 

(3-5) Yqo = Jo^Yqo + K 0j h, 

with Yqo independent of ( Jo,h, Ko^) on the right-hand side of (3.5). 

Remark 3.4. (a) The stationarity condition (3.2) is easy to check. How- 
ever, as the proofs of Theorems 3.1 and 3.3 show, a weaker stationary con- 
dition is the existence of a vector norm || • || and to > such that Jo,t and 
Ko,t satisfy the conditions 

(3.6) Elog||Jo,tJ <0 and Elog+ ||K ,t || < oo. 

By (3.4), .E log || Jo .i || < is equivalent to the requirement that there is a 
strictly positive value of t\ such that the Lyapunov exponent of the i.i.d. 
sequence (J tin , tl („ + i))neNo, tnat is, 

n ! ™o I' Mn-l)Mn ■ ■ ■ J0,tt ||) 

= leNv^ (l0g II J *i(n-l),tin " " " \ti II)) » 
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(which is independent of the specific norm), is strictly negative. As shown by 
Bougerol and Picard [7], provided -Elog + ||«/o,*i|| < °o, Elog + HKo^H < oo 
and a certain irreducibility condition holds, then strict negativity of the 
Lyapunov exponent is not only sufficient, but also necessary for the existence 
of stationary solutions of such random recurrence equations. 

(b) The conditions of Theorem 3.1 imply the conditions (3.6) with the ma- 
trix norm defined as the natural norm r = ||S'- 1 A5|| r , corresponding 
to the vector norm 

(3.7) ||c|| B) r :=||£ -1 c|| rj ceC*. 

Observe, however, that the conditions of Theorem 3.1 are in general not 
necessary for stationarity. For example, using methods similar to those in 
the proofs of Theorems 3.1 and 3.3, it can be shown that for any vector 
norm || • || and for t > 0, 

\\Jo,t\\<\\e Bt \\ + e mit ™p( J2 log(l + (AL s ) 2 ||ea'||)) 

\0<s<t / 

x ||ea'|| ( AL s) 2 - 

0<s<t 

Now if X(B) < 0, then ||e B *|| — > as t — > oo, and (3.6) can be satisfied with- 
out assuming that all the eigenvalues of B are distinct, but choosing ||a|| 
sufficiently small and imposing certain integr ability conditions on L. We 
shall not pursue this argument here because the conditions of Theorem 3.1 
will be sufficient for our purposes. 

The matrices J S) t and the vector K Si j of Theorem 3.3 will be constructed 
explicitly when L is compound-Poisson, and in the general case will be ob- 
tained as the limit of the corresponding quantities for compound-Poisson- 
driven processes. In the compound-Poisson case we shall show that the sta- 
tionary state vector satisfies a distributional fixed point equation which is 
much easier to handle than (3.5). Also, we compare the stationary distri- 
bution of Yqo with the stationary distribution of the state vector when 
sampled at the jump times of the Levy process. This is the content of the 
next theorem. 

Theorem 3.5. (a) Let (Y t )t>o be the state process of a COGARCH(p, q) 
process with parameters B, a and olq. Suppose that the Levy measure vl of 
the driving Levy process L is finite and write the compound-Poisson process 
[L,L]( d ) in the form 

N(t) 

[L,Lp= £ (AL s f=J2 Zi , 

0<s<t i=l 
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where N(t) is the number of jumps of L in the time interval (0,t] and Zi is 
the square of the ith jump size. Let T\ denote the time at which the first jump 
occurs and let Tj, j = 2, 3, . . . , be the time intervals between the (j — l)st and 
jth jumps. Furthermore, let {Tq,Zq) be independent of (Ti,Zj)j e N with the 
same distribution as (Ti, Z\). For i £ No, let 

Ci = (I + Z iea .')e BTi , 
Dj = a Zie 

and T n = Ya=i T% (where Tq := ). Then the discrete time process (Yr n )neN 
satisfies the random recurrence equation 



(3-8) Yr n+1 =C n+1 Yr„+D n+1 , 
Furthermore, for any t> 0, 



n G N . 



,B(t-T N{ t)) 



N(t)-2 

l{AT(t)^0} D Af(t) + C N ( t ) ■ ■ •C A T( t )_jD A r( t )_ i _ 1 

i=0 



(3.9) 



= _B(t-r Ar( f)) 



N(t)-1 



Di + X Ci • • -CjDj+i + C\ ■ ■ ■ C N ^Y 



i=l 



(b) Assume additionally that the conditions of Theorem 3.1 are satisfied. 
Then the infinite sum X)i=o C\ - ■ ■ CjDj+i converges almost surely absolutely 
to a random vector Y, which has the stationary distribution of the sequence 
(Yr n )neN - The stationary state vector Y^ satisfies 



(3.10) 



e BT Y, 



where T is independent of (Tj,Zj)j g pj and has the distribution of T\. Fur- 
thermore, Yqo is the unique solution in distribution of the distributional fixed 
point equation 

(3.11) Yoo^QYoo + R, 

where Yoo is independent o/(Q,R) and 

Q:=e BTo (I + Z ea'), 

H:=a Z e BT °e. 



The fixed point equation (3.11) will play a crucial role in the determina- 
tion of the covariance matrix of Y^, which is studied in the next section. 
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4. Second-order properties of the volatility process. In this section (Yt)t>o 
denotes the state process defined by (2.3), with parameters B, a and ao, 
and driving Levy process L with Levy measure vj_,. The aim of this section 
is to study the autocorrelation function of the volatility process (Vt)t>o- We 
shall write 

/x:= / y 2 dv L (y) and p := / y A dv L (y), 
Jr Jm 

and, if /x < oo (i.e., EL\ < oo), 

(4.1) B := B + /W . 

Observe that -B has the same form as B, but with last row given by (—(3 q + 
ficti, . . . , —(3\ + fia q ). We first give sufficient conditions for the moments of 
Y t to exist. 

Proposition 4.1. Suppose that the eigenvalues of B are distinct, A = 
X(B) < 0, || • || is any vector norm on C q and k G N. Then the following 
results hold. 

(a) IfE\L l \ 2k <oo and E\\Y Q \\ k < oo, 

E||Y t || fe <oo Vt>0. 

(b) If E\L±\ 2k < oo, r G [l,oo], S is a matrix such that S~ l BS is diag- 
onal and 

f ((1 + HS-WSUr y 2 ) k - 1) du L (y) < -Xk, 
Jm. 

then S and r satisfy (3.2) and E\\Y O0 \\ k < oo. In particular, ^(Yoo) exists 
if 

(4.2) EL\<oo and ||5 _1 ea / S'|| r /i < -A, 
and the covariance matrix cov(Yoo) exists if 

(4.3) EL\ <oo and ||S , " 1 ea / ^'|^ p < 2(-A - ||S f - 1 ea / S'|| r .//). 

Furthermore, (4.3) implies (4.2), and (4.2) implies that all the eigenvalues 
of B have strictly negative real parts, in particular that B is invertible and 
fl q / aiix. 

Next, we determine the autocovariance function of the (not necessarily 
stationary) volatility process of Definition 2.1. 
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Theorem 4.2. Let (Vt)t>o be the volatility process specified in Defini- 
tion 2.1, with state process (Y t ) t >o and parameters B, a and oq. Suppose 
that EL\ < oo and that E\\Yt\\ 2 < oo Vt > (as is the case, e.g., if the con- 
ditions of Proposition 4.1 are satisfied). Then, with B defined as in (4.1), 

(4.4) cov(V t+h ,V t ) = a'e Bh cov(Y t ) S L, t,h>0. 

Since we are primarily interested in the stationary volatility process, we 
need to evaluate cov(Y 00 ). However, first we need an expression for E^Y^). 

Lemma 4.3. Suppose that all the eigenvalues of B are distinct and 
that (4.2) holds. Then 

(4.5) E(Y OQ ) = -a »B- l e = ei . 

p q - a\[i 

The following theorem contains the main results of this section. It demon- 
strates that the autocorrelation function of the stationary COGARCH volatil- 
ity process is the same as that of a continuous-time ARMA process. This 
reflects the corresponding discrete-time result that the autocorrelation func- 
tion of a GARCH volatility process is the same as that of a discrete-time 
ARMA process. 

Theorem 4.4. Suppose that the eigenvalues of the matrix B are dis- 
tinct, X(B) < and (4.3) holds. Then the matrix (I ® B) + (B <g> I) + p((ea')(g> 
(ea')) is invertible and the covariance matrix ofY^ is the unique solution 
of 

[(/ ®B) + (B®L) + p((ea!) <g) (ea'))] vec(cov(Y oc )) 

(4.6) 

-OoPqP 



(Pq - fldl 



vec(ee'). 



Let (ipt)t>0 be a stationary CAKMA(q,p — 1) process (as defined in Sec- 
tion 2) with location parameter 0, moving average coefficients a±, . . . ,a p , 
autoregressive coefficients p\ — pa q ,(32 — pa q -\, . . . ,(3 g — a±u, driving Levy 
process L and corresponding state process ((t)t>o- Suppose that E(Li) 2 < oo, 
E(L\) = p and var(Li) = p, and define 



m := pi 8Le Bt ee'e Bt &dt = v&r('i()t 
Jo 



Then < m < 1 and 



cov(Y oc ) = — — rCOv(Co 

(Pq- paiY(l-m) 
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(4.7) 



<*oPip 



(4.8) var(F 0< 

(4.9) E(V V 
(4.10) 



(P g - pai) 2 (l - m) Jo 

a oPg 171 

{{3 q — fiai) 2 1 — m ' 
j3 q - \ia\ ' 



re Bt ee'e B,t dt, 
Jo 



Pq ~ potl 

If (Vt)t>o is the stationary COGARCH volatility process, then 

OtnPn 

(4.11) cov (V t+h ,V t ) = —- 0H J cov(^ +fe ,^), t,h>0, 

\jjq /JiCtl) yi m) 

showing, in particular, that V has the same autocorrelation function as ip. 
If the eigenvalues Ai,...,A 9 of B are also distinct, and a(z) and b(z) are 
the characteristic polynomials associated with a and B, then 

cov(V t+h ,V t ) - 



(P q - pai) 2 (l - m) 
(4-12) 

xVfttMe^, t,h>0, 
where b' denotes the derivative of b. 



5. Positivity conditions for the volatility. For the definition of the 
COGARCH price process dGt = yfVtdt to make sense, it is necessary that 
Vt be nonnegative for all t > 0. The following theorem gives necessary and 
sufficient conditions for this to occur with probability 1. 

Theorem 5.1. (a) Let (Y t )t>o be the state vector of a COGARCH(p, q) 
volatility process (Vt)t>o with parameters B, a and ao > 0. Let 7 > — cxq be 
a real constant. Suppose that the following two conditions hold: 

(5.1) a'e B *e>0 Vi>0, 

(5.2) a'e Bi Y >7 a.s. Vt>0. 

Then for any driving Levy process, with probability 1, 

(5.3) V t >a o + -y>0 Vt>0. 
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Conversely, if either (5.2) fails or (5.2) holds with 7 > —ao an d (5-1) fails, 
then there exists a driving compound- Poisson process L and to>0 such that 
P(V t0 < 0) > 0. 

(b) Suppose that all the eigenvalues of B are distinct and that (3.2) 
and (5.1) both hold. Then with probability 1 the stationary COGARCH(p, q) 
volatility process (Vt)t>o satisfies 

V t >a >0 Vf>0. 

For the stationary COGARCH volatility process or for the process with 
Yo = 0, the condition (5.1) alone is sufficient for almost sure nonnegativity. 
The expression a!e Bt e is in fact the kernel of a CARMA process with au- 
toregressive coefficients bi,...,b q and moving average coefficients ai,...,a q . 
Results that pertain to nonnegativity of a CARMA kernel were recently ob- 
tained by Tsai and Chan [28]. We state their results in the next theorem in 
the context of COGARCH rather than CARMA processes. Statement (e) 
below was also obtained by Todorov and Tauchen [27]. Recall that a func- 
tion (j> on (0,oo) is called completely monotone if it possesses derivatives of 
all orders and satisfies (-l) n {d n <j)/dt n ){t) > for all t > and all n £ N . 

Theorem 5.2. Let B and a be the parameters of a COGARCH(p, q) 
process. If X(B) < 0, and ct\ > 0, we have the following results. 

(a) For the COGARCH(p, q) process, (5.1) holds if and only if the ratio 
of the characteristic polynomials a(-)/b(-) is completely monotone on (0,oo). 

(b) A sufficient condition for (5.1) to hold for the COGARCH(l,g) pro- 
cess is that either (i) all eigenvalues of B are real and negative or (ii) if 
(Aj 15 Aj 1+ i), . . . , (Aj r , Aj r+ i) is a partition of the set of all pairs of complex 
conjugate eigenvalues of B (counted with multiplicity), then there exists an 
injective mapping u : {1, . . . ,r} — > {1, . . . , q} such that X u tj\ is a real eigen- 
value of B satisfying A u (j) > ^(A^). 

(c) A necessary condition for (5.1) to hold for the COGARCH(l,g) pro- 
cess is that there exists a real eigenvalue of B not smaller than the real part 
of all other eigenvalues of B. 

(d) Suppose 2 < p < q, that all eigenvalues of B are negative and ordered 
as in Definition 2.3, and that the roots jj of a(z) = are negative and 
ordered such that 7 P _i < • ■ • < 71 < 0. Then a sufficient condition for (5.1) 
to hold for the COGARCH(p, q) process is that 

k k 

£7*<E A * Vfc€{l,...,p-1}. 

i=l i=l 

(e) A necessary and sufficient condition for (5.1) in the COGARCH(2, 2) 
case is that both eigenvalues of B are real, that cti > and that a\ > 
-a 2 X(B). 
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Although characterization (a) may be difficult to check in general, it gives 
a method for constructing further pairs (a, B) for which (5.1) holds, since 
the product of two completely monotone functions is again completely mono- 
tone. 



6. The autocorrelation of the squared increments. In Section 4 we inves- 
tigated the behavior of the autocorrelation function of the volatility process. 
Since one of the striking features of observed financial time series is that the 
returns have negligible correlation while the squared returns are significantly 
correlated, we now turn to the second-order properties of the increments of 
the COGARCH process itself. We therefore assume that V is strictly sta- 
tionary and nonnegative, and define, for r > 0, 

G^ := G t+r - G t = f Vv s dL s , t>0. 

J(t,t+r] 

It is easy to see that (G[ )t>o is a stationary process. Let /i and B be defined 
as in Section 4. We then have the following theorem. 



Theorem 6.1. Let B, a and ao be the parameters of a COGARCH(p, q) 
process whose driving Levy process has mean zero. Suppose that the eigenval- 
ues of B are distinct, that (4.2) and (5.1) hold, and that V is the stationary 
volatility process. Then for any t>0 and h>r>0, 

(6.1) E(G^) = 0, 

(6.2) E{{G?f) = *&-E(Ll), 

p q - m\ 

(6.3) cov(Gf ) ,G^) = 0. 
If in addition (4.3) holds, then 

(6.4) cov((G| r) ) 2 , (G[%) 2 ) = a'e^H r , h > r, 
where 

H r := E{L\)B-\I - e~ Sr ) cov(Y r , G 2 r ). 



The autocovariance function (6.4), like that of the CARMA process with 
parameters B and a, is a linear combination of terms of the form e Xjh , 
j = 1, . . . , q, where . . ,\ q are the eigenvalues of B. 



18 



P. BROCKWELL, E. CHADRAA AND A. LINDNER 



7. An example. In this section we illustrate the properties established 
above using the COGARCH(l,3) process driven by a compound-Poisson 
process with jump rate 2 and normally distributed jumps with mean and 
variance 0.74. The COGARCH coefficients are ao — a i = 1, /3i = 1-2,(32 — 
0.48 + 7r 2 and (3% = 0.064 + 0.47T 2 , from which we find that the eigenval- 
ues of B are —0.4, —0.4 + ttz and —0.4 — ixi. With S defined as in (3.1), 
||S ,_1 ea'5||2 = 0.21493 and it is easy to check from this that the conditions 
(4.2) and (4.3) are satisfied. Condition (b)(ii) of Theorem 5.2 also implies 
that the volatility process is nonnegative. 

The eigenvalues of the matrix B = B + /W are -0.25038, -0.47481 + 
3.14426* and -0.47481- 3. 14426z. From (4.12) we conclude that the autocor- 
relation of the volatility in this case is a linear combination of exp(— 0.25038t), 
and a damped sinusoid with period approximately equal to 2 and damping 
factor exp(-0.4748U). 

The top graph in Figure 1 shows the values at integer times 101, . . . , 8100 
of a simulated series (Gt) with the parameters specified above, Yo = 
(1, 1, 1)' and G(0) = 0. The second graph shows the differenced series (Gt+i — 
G ! t)t=ioo,...,8099 and the last graph shows the volatility (<r i 2 )t=ioi,...,8ioo- 

As is the case for a discrete-time GARCH process, the increments (Gt+i — 
Gt) exhibit no significant correlation, but the squared increments ((Gt+i — 
Gt) 2 ) have highly significant correlations as shown in the second graph of 
Figure 2. The first graph in Figure 2 shows the sample autocorrelation func- 
tion of the volatility process at integer lags. This too is highly significant 
for large lags, reflecting the long-memory property frequently observed in 
financial time series. As expected from the remarks in the first paragraph 
above, it has the form of an exponentially decaying term plus a small damped 
sinusoidal term with period approximately equal to 2. 

8. Proofs for Section 3. We start by proving Theorem 3.5, since (3.9) 
will be needed in the proof of Theorems 3.1 and 3.3. 

Proof of Theorem 3.5. (a) It follows from (2.3) that Y t satisfies 
dYt = B~Yt dt for t £ [T n , r n +i), so that 

(8.1) Y^e^'-^Yrv, t e [F n ,T n+1 ),n G N . 

At time T n+ i a jump of size e(«o + a'Yr„ +1 _)Z n+ i occurs, so that 

Y r „ +1 = Yr n+1 _ + e(a + a'Y r „ +1 -)Z n+1 
= (I + Z n+ iea')Yr n+1 - + a Z n+ ie 
= C n+ iYr n +D n+ i, n G No, 
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Fig. 1. The simulated compound-Poisson driven COGARCH(l,3) process with jump 
rate 2, normally distributed jumps with mean and variance 0.74, and coefficients 
oto = on = 1, /3i = 1.2, 02 = 0.48 + 7r 2 and /3s = 0.064 + 0.4-7T 2 . The graphs show the process 
(Gt) sampled at integer times (top), the corresponding increments ((Gt+i — Gt)) (center) 
and the corresponding volatility sequence (Vt = of ) (bottom). 
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Fig. 2. The sample autocorrelation functions of the volatilities (Vt) (left) and of the 
squared COG ARCH increments {{Gt+i — Gt) 2 ) (right) of a realization of length 1,000,000 
of the COGARCH process with parameters as specified in Figure 1. 



which is (3.8). Solving this recursion gives 

n-2 

Y r „ = D n + Cn ■ ■ ■ C re -iD re _i_i + C n ■ ■ ■ dY , n G N, 

8=0 

and the first equality in (3.9) follows from this and Y t = e B ^~ rjv <''^Yr JV(t) . 
The second equality in (3.9) is a consequence of the fact that the infinite 
random element (N(t),T^u\,C^(t)^N(t)) ■ ■ ■ >CijDi,0,0, . . .) has the same 
distribution as (N(t),T N ^,Ci,T)i, . . . , CWt)> T) N ^, 0, 0, . . .); indeed, for any 
n G No and c > 0, the random vectors (Ci, Di), . . . , (C n , D n ) are i.i.d. and 
depend on the restriction {N(t) = n, ^N(t) > c } only in terms of J2?=i^i 
and T n+ i, but not on the Tj, i = 1, . . . ,ra, individually. 

(b) Let S be such that S~ 1 BS =: A is diagonal and define the vector norm 
|| c lls,r — H5 c|| r as in (3.7), so that the associated natural matrix norm is 
\\A\\ B , r = WS^ASWr. Then we have for t > 0, 

(8.2) ||e Bt || B , r = ||Se A *5- 1 || B)r = ||e At || r = e M . 

This gives ||Ci||b )7 . < (1 + Zi ||ea'||s,r)e ATl and ||Di||_B )r = CKo||e||s,r Zi, so 
that, using Vu,n([x, oo)) = vl{v G M. : \y\ > y/x} for x > 0, 

^log||Ci||B,r < AE(Ti) + £log(l + Z 1 ||ea , || Bir ) 

= — 7^ + —m\ I log(l + ||ea / || Bir y 2 )^ L (y)<0 
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by (3.2) and 

Elog + (Z 1 ) = ^— f log+(y 2 )du L (y)<oo. 

From the general theory of random recurrence equations, this implies the 
almost sure absolute convergence of J2i=o C\ - ■ ■ CjDj + i to Y, which has the 
stationary distribution of (Yr n ) n eN (see, e.g., [7]). 
To prove (3.10), for m 6 N, let 

m—l 

Y m := 2J Ci • • ■ CiDi + \ + C\ ■ ■ ■ C m Yo 

and 

Y t , m :=e B ^- TN (^Y m , t>0. 

Since the random variable (t — r^vft)) is asymptotically independent of T\, 
Zi, . . . ,T m , Z m (for t — > oo, m fixed), it follows that e B<yt ~ TN ^ is asymp- 
totically independent of Y m and, hence, Y^ m converges in distribution to 
e BT ~Y m as t — > oo, where T is exponentially distributed with parameter 
vl(E) (e.g., [26], Section 7.4.4) and is independent oiT\,Z\,..., T m , Z m and, 
hence, can be chosen to be independent of (Tj)j e pj, (^)ieN ( as in the state- 
ment of the theorem). Moreover, e BT Y m converges almost surely, hence in 
distribution to e BT Y, as m — ► oo. Denote by Y t the expression in the lower 
line of (3.9). Then (3.10) and, in particular, the existence of the limit vari- 
able Yqo in the compound-Poisson case follow from (3.9) and a variant of 
Slutsky's theorem (e.g., [11], Proposition 6.3.9), provided 




Since ||e B (*- r *w)|| Bjf .<l and l {JV(t) ^ 0} Di + e£i 1 Ci • • • dD i+1 + d • • • 
C^v( t )Yo — Y m converges almost surely, hence in probability as t — > oo to 
J2i^=m Ci " • " CjDj + i — (7 X . . . Cm Yo, which itself converges almost surely to 
asm-> oo, (8.3) is true and (3.10) follows. That Yoo satisfies (3.11) is clear 
from (3.10); that it is the unique solution follows from Elog \ \Q\\B,r < and 
Elog + ||R||s, r < oo. □ 

The proof of Theorem 3.5(b) already showed the existence of the limit 
variable Y m for the case of a driving compound-Poisson process. Neverthe- 
less, this existence will be reestablished in the proof of Theorems 3.1 and 3.3 
for the general case, making use of Theorem 3.5(a) only. We shall use an 
approximation argument and introduce the following notation: 
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Definition 8.1. Let L be a Levy process. Then for any e > 0, the y/e- 
cut Levy process (if )oo is defined by 

L [ t ] := ALs ' 
0<s<t,|AL s |> v / i 

If (Y t )t>o is a state process of a COGARCH(p, g) process driven by L, then 
the state process of the COGARCH(p, q) process with the same parameters 
and starting vector but driving Levy process (i^ )t>o will be denoted by 

The quadratic covariation of is given by 

[lW,lW] 4 = [lW,lW]« 

0<s<t,|AL s | 2 >e 

In particular, the corresponding COGARCH volatility is driven by a compound- 
Poisson process. With this notation, we have the following lemma: 

Lemma 8.2. Let (Y t )t>o be the state process of a COGARCH(p, q) pro- 
cess. Then Y^ converges in ucp to Yj as e — > 0. 

Proof. This is an easy consequence of perturbation results in stochastic 
differential equations: recalling the definition of prelocal convergence in H p , 

1 < p < oo, as in [24], page 260, it is easy to see that [L( e \L^} converges 
prelocally to [L, L] w in H p , 1 < p< oo, ase^O (e.g., with stopping times 
Tfc = k). The claim then follows from Theorems 14 and 15 of Chapter V in 
[24]. □ 

Proof of Theorems 3.1 and 3.3. We shall first concentrate on (3.3) 
and (3.4), and then prove Theorem 3.1 and the rest of Theorem 3.3 simul- 
taneously. Let e > and assume the representation 

N e (t) 
i=i 

where L^ is the y/e-cut Levy process of Definition 8.1. Define cf ) and Dj e) 
similarly as in Theorem 3.5. Furthermore, let 

j(e)._ B(*-r« ) Me) ___ r ( £ ) 



Kg:= e B(t - r ^ W ) 



N e (t) 
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N e (t)-2 

, \ " M £ } r 1 ^ n( £ ) 

+ U N e (t)'" U N E <t)-i U 



N s (t) N E {t)-i N £ (t)-i-l 

i=0 



Then, by Theorem 3.5(a), 



(8.4) YW = J$Y + Kg. 

Prom the previous lemma we know that Y^ converges in ucp to Y^ as 
e — > 0. Since this is true for any starting value Yq, it holds in particular for 
Yo = 0, and from (8.4) it follows that K e ] converges in ucp to some Ko,t as 
e — > 0. Hence, again from (8.4), it follows that for arbitrary Yo, 

J o,t Y o— Y t — tt-o.t Y 4 ~~ ^-o.t ase^U. 

Since this holds for arbitrary Yo, we conclude that JqJ converges in ucp to 
some Jo,t as e — > 0. From (8.4) it then follows that 

Yt = Jo,*Yo + Ko,t- 

By starting at an arbitrary time s instead of at time 0, we obtain (3.3). For 
example, j\ £ } is given by 

j(e) = e B ^l {t) ) CNe{t) • • • C Ne{s)+2 {I + Z Ns{s)+1 e a ')e B{T ^m- s \ 

0<s<t, 

giving (3.4). The independence and stationarity assertions on (J Sj j,K S] t) are 
clear, since J Sj t and K Sji are constructed only from the segment {L u ) s<u <_ t 
of the Levy process L. 

Now assume that all eigenvalues of B are distinct and that (3.2) holds. 
Applying (8.2) to J ° 

H7 {e) ll < ||e B{t - r ^w)|| \\C {£) II •■•ll<7^|| 
ll^O.i \\B,r — H e II B,r II u jV e (t) II -B,r W^l \\B,t 



Ne(t) 

<e A(t - r -e( t )) n ((l + Z^||ea'|| B)r ) e ^ e) - r a)) 
i=i 



.5) = eA * ex p(^E loga+zfis-wsiyj 

.6) <e At exp( ^ log(l + (AL s ) 2 ||5~ 1 ea / 5|| r )J. 

\0<s<i / 
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Since || Jo,t\\B,r < limsup £ ^ ||Jot we conclude that 

(8.7) log||J ,t||B,r<At+ £ log(l + (AX,) 2 !!^-^^^!^), 

0<s<t 

giving 

£log || Jo,t\\B,r <t(\ + jf log(l + WSea'S^Wy 2 ) du L {y)^j < 

by (3.2) (see, e.g., [24], Chapter I, Theorems 36 and 38). This is the left-hand 
inequality of (3.6). To show that Elog + ||Ko,t||s,r < oo, observe that 



IK (£) II <p 

l-^O.t \\B,r — e 



( r %wh {Ne(t)m a \\e\\ Bjr Z$ 



(t) 



N e {t)-2 



(8.8) 



(0 



+ ao||e||B, r £ e^" 1 ^^)— ^l + Z^Hea'H^) 



x--.x(l + 4 £ ) .||ea'|| B , r )Z (£; 



N e (t)-i-l 



< ao\\e\\B,rl{N e (t)^0}ZN E (t) 
N e (t)-2 

+ ao||e||B,r £ ex P 

i=0 lQ<s<t 



£ log(l + (AL s ) 2 ||ea'|| B , 



JV £ (t)-i-l 



< ao||<S 1 e|| r exp 



(8.9) 



£ log(l + (AL a ) 2 ||5-Ws r || J 

0<s<t 



x £ (AL S ) 2 . 

0<s<t 

From this it follows that 

log ||Ko,t||B,r < log(a ||'S ,_1 e||r) 

+ E log(l + (AL s ) 2 ||5^ 1 ea / > S|| r ) + log[L,L]| d) . 

0<s<t 

The expectation of the second summand is finite as shown above and -E(log[L, 
L]i) <°° since I(i,oo) lo Sxdu [LiL] (x) =/ R \ [ _ 1 ^logx 2 diy L (x) < oo, showing 
the right-hand inequality of (3.6). 

Let (J n ,K n ) n6 N be an i.i.d. sequence with distribution (Jo^Ko,].), inde- 
pendent of L and Yo Let 7 € [0, 1) and n £ N. Then it follows from (3.3) 
that 

n-2 

Yyi-|_-y K n -|_>y_l )n -)_>y + ^ ' J n +^— \ ,71+7 ' ' ' ^Tl+J— % — 1 ,71+7" j K n _|_-y _ j _ 2 ,71+7 — 1— 1 

t = 



CONTINUOUS-TIME GARCH PROCESSES 



25 



"Ai-Hy— 1,71+7 ' ' ' ^y^+l^T 1 

d 

— Ki + J\ ■ ■ ■ JiKj + i + Ji • • • J n Y 7 
1=1 

= :G n + iJ n Y 7 , say. 

Since -Elog || Ji||s,r < an d -Elog" 1 " ||Ki||b )7 . < °o, it follows from the general 
theory of random recurrence equations (e.g., [7]) that H n converges almost 
surely to as n — > oo and that G n converges almost surely absolutely to 
some random vector G as n — ► oo. Since Y has cadlag paths, it follows that 
sup 7g [ 0) i) ||Y 7 ||b |7 . is almost surely finite. Hence 

lim sup ||i? n Y 7 ||s r = a.s., 

76[0,1) 

and it follows that Yt converges in distribution to Yoo := G as t — > oo. That 
Yqo satisfies (3.5) and is the unique solution is clear by the theory of random 

recurrence equations. Equations (3.5) and (3.3) then imply that if Yo = Y m , 

then Yt = Y m for all t > 0, showing strict stationarity of (Y t )t>o since it is 
a Markov process. □ 

9. Proofs for Section 4. To prove Proposition 4.1, we will show that 
the state process (Y t ) t >o can be majorized by the state process of a 
COGARCH(l,l) process, for which we can apply the moment conditions 
of Kliippelberg, Lindner and Mailer [19]. We further show that under the 
conditions of Theorem 3.1, the stationary distribution Y m can be approxi- 
mated by stationary distributions of compound Poisson-driven COGARCH 
processes and that there is a majorant for this approximation. This will allow 
us to restrict attention to compound Poisson-driven processes when calculat- 
ing autocorrelations, the general case following from Lebesgue's dominated 
convergence theorem. This is the content of the next lemma. 

Lemma 9.1. Let (Y t )t>o be the state process of a COGARCH(p, q) pro- 
cess with parameters B, a and cto > such that all eigenvalues of B are 
distinct and that A = X(B) < 0. Let r € [1, oo] and let S be such that S~ 1 BS 
is diagonal. Denote by || • ||_b r the vector norm defined in (3.7). Furthermore, 
denote by (Yt)t>o the state process of a COGARCH(l,l) process satisfying 
(2.3) with the parameters (B,a,ao) replaced by (A, ||ea'||.B,r> Q!o||e||B,r) an d 
initial state vector Yq := ||Yo||s r . Then 

(9.1) ||Y t || B| r<Y't, t>0. 

If (3.2) is satisfied for this r, then there exist versions ofY^ and Y M such 
that 

(9.2) HYooIIb^Yoo. 
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Furthermore, if (Y^ )t>o is the process defined in Definition 8.1 for e > 0, 
then versions of Y^j can be chosen such that ||Y^'||s r < Y m for all e > 
and Yoo, as e — > 0. 

Proof. We use the notation and setup of the proof of Theorems 3.1 
and 3.3. Let e > and define a COGARCH(l, 1) state process Y similarly 
as above (with respect to Y^). Let Jq\ and Kg] be defined similarly as 
J^fl and Kq £ ] (with respect to Y ). Then it is easy to see that JqJ and 
Kq £ ] are the right-hand sides of (8.5) and (8.8), respectively. In particular, 

\\Jot \\B,r < Jq 6 } an d W^o} \\s,r < Kq^, and since Jq\ and K ^ converge in 
ucp as e — > to some J 0j j and K .t such that 

Yt = J o,tYo + Ko,t, 

it follows that ||Y t ||s,r < Y t for fixed t > 0, giving (9.1). 

Similar quantities such as and J Sj j can be defined when going from 
time s to time t, and similar results hold. Let := ao|| e ||s r + || ea/ ||s r Y^_ 
be the COGARCH(l, 1) volatility corresponding to Y (e) . Define 



X t :=-Xt- logCl + CAL.fllea'llB,, 



0<s<t 



X\ £) :=-\t- log(l + (AL s ) 2 ||ea'|| Bir ). 

0<s<t,(AL a ) 2 >e 

Then it follows from Theorem 2.2 and (1.4), that 

= (Vo - a \\e\\ B , r \ J\ x ° £) ds)e- x t S \ 

Thus we have Jq \ = e~ x t ' and obtain another formula for K Q e ] , namely 

_v-00 , f f _fx (£) -X (e, l , 

ao|| e l|sre * — ao|| e l|BrA / e 1 * s ; as — cto|| e lls 

JO 

From this it can be seen that \ and K Q e ] are bounded by Jqj = e - ^* and 



t7^ £ ) ll 'li— 1 
K ,t = ll ea II R r 



K ,t = ||ea'|| B ^.Q!o||e||s )7 . 
respectively. Now define the versions 

oo 

Yqo := ^2 ^0,1 " Ji-l,iK-i,i+l) 



e~ Xt - A 



Jo 



i=0 
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oo 

V(<0 — V 7 (e) ••• T (e) K (e) 

i=0 
oo 

Yqo := ^ Jo,i • • • Ji-i^Kj^+i. 

i=0 

In the proof of Theorems 3.1 and 3.3 we have seen that (3.2) implies that 
the sum defining Yoo converges almost surely. This then gives the claim, 
since 

II T II II II <r 7 
\\<Ji—l,i\\B,ri II u i—l,i\\B,r — i— 

IIK^+iIIb^, ||K^^ +1 || Bjr < Ki )i+ i, 

and J>f.\ i and K^ + i converge in probability to J%-\-i and Kj^+i as e — > 0, 
respectively. □ 

Proof of Proposition 4.1. All assertions apart from the implication 
"(4.2) ==^ X(B) < 0" follow immediately from Lemma 9.1 (observing that 
the existence of is independent of the specific matrix norm) and 

the corresponding properties of the COGARCH(l,l) process; see Section 4 
in [19]. That (4.2) implies X(B) < is a consequence of the Bauer-Fike 
perturbation result on eigenvalues, stating that for every eigenvalue Xj of B 
we have 

min |A;-L| < \\S- 1 (B-B)S\\ r = ii\\S- 1 ea!S\\ r 

i=l,...,q 

(see, e.g., Theorem 7.2.2 and its proof in [17]). □ 

Proof of Theorem 4.2. Since for fixed t, almost surely V t = V t+ = 
ao + a' Yt , we obtain 

(9.3) cov(V t+h , V t ) = a' cov(Y t+h , Y t ) a. 

For ease of notation, we will assume that t = 0. Let := Jq^ and K/j := K 0i h 
as constructed in the proof of Theorem 3.3. Then, using that \\e Bt \\ < e^H* 
for any vector norm || • ||, it follows as in the proof of (8.6) that 

(9.4) £?||J fc ||<ell B ll*^exp| V log(l + (AL s ) 2 ||ea'||) ) } < oo 



;^exp Yl 

I \0<s<h 



by Kliippelberg, Lindner and Mailer ([19], Lemma 4.1(a)). Using that 
J/ t Yo + Kfc, we conclude that S||K^|| < oo and that 

E(Y h Y' ) = E(E(Y h Y' \J h ,K h )) 

= E(J h E(Y Y' ) + K h E(Y' Q )) 

= E(J h )E(Y Y' ) + E(K h )E(Y' ). 
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On the other hand, 

E(Y h ) E(Y' ) = E(J h )E(Y )E(Y' ) + E(K h )E(Y' ), 

so that cov(Y/ l ,Yo) = E{ J/J cov(Yo) and (4.4) will follow from (9.3) once 
we have shown that 

(9.5) E(J t ) = e Bt , t>0. 

To do that, it suffices to assume that [X, £]* is a compound- Poisson process. 

The general case then follows from the fact that as defined in the proof of 
Theorem 3.1 converges to Jt in L 1 as e — > 0, since it converges stochastically 
and since there is an integrable majorant by (9.4) and its proof. So suppose 

that [L, L]t = Si=i %i is compound-Poisson with intensity c > and let 
d = (I + Z i ea / )e B ( r *- r <- 1 ). Then, for < s,t, it follows from (3.4) and the 
independence of Jo )S and J s , s +t that 

E{J s+t ) = E(J s )E(J t ). 

It is easy to see that E(J t ) is a continuous function in t € [0,oo). Further- 
more, E{Jq) = I and we conclude that (E(J t ))t>o is a semigroup. We shall 
show that its generator A j satisfies 

(9.6) Aj:=Yiui-{E{J t )-I) = B+ f y 2 du L {y) ea' = B . 

This then implies (9.5), since E(J t ) = e tAj (see, e.g., [16], Proposition 2.5). 
To show (9.6), write 

Jt-e l {N (t)=o}+e ^il {JV(t)=1} 

+ e B( - t - T ^ ) C N{t yC 1 l {Nm2} . 

Since N(t) is Poisson distributed with parameter ct, we have P(N(t) = k) = 
e- ct (ct) k /(k\). Then by (9.4), 

E{e B ^- T ^C my --C 1 l {m > 2) ) 

^elW^^exp^loga + Zllea'll^l^)^ 

= e^^exp^ log(l + Z||ea'||)^ |iV(t) > 2^ P{N{t) > 2) 
/N(t)+2 



(9.8) <e mt E 



^expj^ J2 log(l + Z l ||ea / ||)JJp(iV(t)>2) 



elWECCl + ZiHea'IDCl + Zallea'l 



CONTINUOUS-TIME GARCH PROCESSES 29 
x i^exp^ J2 + (AL s ) 2 ||ea / ||)^P(iV(t) > 2) 

= o(t) as t — ► 0, 

since P(N(t) > 2) = o(i) as t — > 0. Furthermore, since Ti is uniformly dis- 
tributed on (0,t), conditional on N(t) = 1, it follows that 

S(e s (*- ri )C 1 l {JVW=1} ) 

= E{e B{t ~ ri \l + Ziea')e Bri |iV(i) = l)P(N(t) = 1) 

= /* e B ^ (I + ^(ZOeaOe 55 — e~ ct ct. 
Jo t 

Since sup 0<s<4 \\e Bs — I\\ converges to as t — > 0, we conclude that 

liml£J( e s (*- r i)c 1 l {JV(t)=1} ) = (J + £(Zi)ea')c. 
Now (9.7) and (9.8) give (9.6), since 

Jim ^ (Jt) ~ 7 = lim eBte " Ct ~ 1 + c(I + E(Z 1 )ea!) 

= -d + B + c(I + E(Z 1 )ea / ) = B. □ 

We now need the following lemma: 

Lemma 9.2. Let T be exponentially distributed with parameter c, and 
suppose that \{B) < 0. Let 

M :=E{e BT ®e BT ). 

Then 

(9.9) E(e BT ) = (I-c- 1 B)-\ 

(9.10) M~ l = L e -{I® (c _1 5)) - ((c~ l B) ® I). 

Furthermore, (I <g> B) + (B<S) I) is invertible and, for any real (q x q) matrix 
U, the unique solution of ((I (8) B) + (B <8> J)) x = vec(C7) is given by 

(9.11) x = vecf- f°° e Bt Ue B,t dt). 



o 



Here, we denote by I the (q x q) identity matrix and denote by I q 2 the 
(q 2 x q 2 ) identity matrix. 
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Proof. Equations (9.9) and (9.10) follow by simple calculations and a 
diagonalization argument, while invertibility of (I ® B) + {B®I) and (9.11) 
is a consequence of Lyapunov's theorem for the solution of Lyapunov equa- 
tions (see, e.g., Section 9.3 in [15]). □ 

Proof of Lemma 4.3. Suppose first that the Levy measure of L is 
finite, and let Q and R be as in Theorem 3.5(b) [writing (T, Z) instead of 
(Tq,Zq)]. Then, by Lemma 9.2, 

E{Q) = (I-c~ l By\l + E(Z)ea.'), 

E(K) = a E(Z)(I - c~ 1 By 1 e, 

and (3.11) gives 

(I-E(Q))E(Y 00 )=E(K). 

Furthermore, 

(I - c" 1 J B)(/ - E(Q)) = [(I - c- x B) -I- E(Z)ea'} 

= --(B + /iea ; ), 

c 

giving 

£(Yoo) = -c(B + fiea')- 1 ^ - c~ 1 B)E(R) = -a n(B + //ea') _1 e. 

Denoting u = (u±, . . . , u q )' := (B + /iea') _1 e, it is easy to see that « 2 = • • • = 
u q = and u\ = l/(«i/i — f3 q ). In the case when vl is infinite, the result 
follows from Lemma 9.1, using that 7^ is an integrable majorant by (4.2). 

□ 

Proof of Theorem 4.4. By Lemma 9.1 and the dominated conver- 
gence theorem, for showing (4.6) it is sufficient to assume that [L, L] is 
a compound-Poisson process. Hence, let Q and R be as in Theorem 3.5, 
writing (T,Z) instead of (Tq,Zq), where T is exponentially distributed with 
parameter c > 0. Then 

^(YooY^) - EiQY^Y'^Q') 

(9.12) 

= ^(QYooR) + E(KY' 00 Q') + E(RB!) 
by (3.11) and all these expectations exist by (4.3). Now 
EiQY^Y'^Q') = EiEiQY^Y'^Q'lQ}) 

= E(E[QE{Y 00 Y' 00 )Q'\T]) 

= E(e BT E[(I + Zea / )i?(Y 00 Y^)(J + Zae')]e B ' T ). 
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Using that vec(AiA2A 3 ) = (A' 3 <g> A{) vec(A 2 ) for matrices Ai,A 2 and A3, it 
follows with M as in Lemma 9.2 that 

vec^QYooY^Q')) 

= Mvec(£((I + Zea!)E{Y O0 Y' oo )(I + Zee'))) 

= M{E((I + Zeal) + Zea'))) vec(^(Y oc Y^)) 

= M(I q2 + E(Z)((ea') ® I) + E(Z)(I <g> (ea')) + E(Z 2 )((ea') ® (ea'))) 

x vec^Y^)). 

Similar expressions can be obtained for vec(£'((5Y 00 R / )), vec(i?(RY^ (5 / )) 
and vec(-E'(RR / )), and we obtain from (9.12) that 

[I q2 - M(I q 2 + E(Z)((ea!) ® /) + (ea')) + £(Z 2 )((ea') ® (ea')))] 

x vec(i?(Y 00 Y^)) 

= Mvec[a§£(Z 2 )ee' + a (E(Z)I + J B(Z 2 )ea')^(Y 00 )e / 

+ a e£(Y;j(£(Z)/ + £(Z 2 )ae')]. 

Multiplying this equation by cAf , using (9.10) and (4.5) as well as p = 
cE(Z) and p = cE(Z 2 ), we obtain 

- [(/ g> (B + pea')) + ((-B + ilea!) <8> J) + p((ea') <g> (ea'))] vec( J E(Y 00 Y^)) 
= vec[aopee' — a (/j,I + pea!)n(B + pea / ) _1 ee / 
— ao ee '(-B' + pae / ) _1 p(pl + pae')] . 

Adding to this 

[(/ ® B) + (B ® I) + p((ea') (ea'))] vec(£;(Y 00 )£;(Y^)) 

= vecfS EiY^EiY'^) + £(Y 0O )£(Y^)B' + pea'^Y^i^Y^Jae'] 
= a 2 , vec^ee'iB')- 1 + fj?B~ x es! + ppWi^W^TW] 
on both sides results in 

-[(/ ® B) + (B ® /) + p((ea') (ea'))] vec(cov(Y )) 
= agp[l - p(a / 5" 1 e)] 2 vec(ee') 

= 7~R \2 VeC ( ee h 

which is (4.6), where we used (4.5) in the last equation. 

Now let A := (I ® B) + (B <g> I) and x := vec(cov(Yoo)). By Proposition 4.1 
and Lemma 9.2, A is invertible. Observe that the matrix p((ea') ® (ea')) has 
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nonzero entries only in the last row. Denote this row by c'. Furthermore, set 
7 := pa>Q(3q([iat — Pq)~ 2 ■ Then (4.6) can be written as 

Ax. + (c'x)e <? 2 = — 7e„2. 

We know already that a solution to this equation exists. Suppose there are 
two solutions, call them xi and X2- Then Ax\ = — (7 + c'xi)e 9 2 and Ax 2 = 
— (7 + c'x2)e (? 2. Denoting the unique solution of Ay = —ne q 2 by y(n), n £ M, 
it follows that xi = y(7 + c'xi) and X2 = y(7 + c'x2). Since xi 7^ 7^ X2, 
this implies 7 + c'xi 7^ 7^ 7 + c ; X2, and using the linearity of the solution 
y(n) in n, it follows that there is k 7^ such that X2 = rexi. Thus we have 
Ax\ = — (7 + c'xi)e (? 2 and kAx\ = — (7 + Kc'xi)e (? 2, and this is only possible 
if k = 1, so xi = X2. So the solution of (4.6) is unique, implying that the 
matrix A + p((ea') (g> (ea')) is invertible. 

By (9.11), the solution y(n) of Ay = —ne q 2 is given by 

(9.13) y(n) = vec(n J™ e Bt ee' e B ' 1 dt^j . 
This gives 

POO ~ ~ 

cov(Yoo) = (7 + c'vec(cov(Y 00 ))) / e Bt ee'e Bt dt. 

Jo 

Since both cov(Y DO ) and f£° e Bt ee'e B/t dt are positive semidefinite, it follows 
that 7 + c' vec(cov(Yoo)) > 0. By Brockwell [9], the stationary CARMA state 
vector has covariance matrix 

[•00 ~ ~ t 

cov(Coo) = p / e Bt ee'e B * dt, 
Jo 

so that there is u > such that 

(9.14) cov(Yoo) = ucov(Coo)- 
Inserting (9.14) into (4.6) and using (9.13) shows 

-up vec(ee') + up 2 vec ^ea' J e Bt ee'e B,t dt ae'^j 
-aoPqP 



vec(ee'), 



so that 



Since u > and ao,f3 q j^ 0, it follows that < m < 1 and that 



-u{l — m) vec(ee ) = — ^—^ vec(ee ) 



*lP 2 q 



(P q - pai) 2 {l-mY 
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giving (4.7). This implies (4.8), using = ao + a 'Yoo, and (4.9) follows 
from (4.5). Finally, 

/•oo — roo — 

E(ip 00 )=a'E e Bt edL t = n &'e Bt edt = -fi&'B^e, 
Jo Jo 

giving (4.10), and (4.11) and (4.12) are direct consequences of (4.4), (4.7) 
and the autocovariance function of a CARMA process (see [9]). □ 

10. Proofs for Section 5. 

Proof of Theorem 5.1. (a) Suppose that (5.1) and (5.2) both hold. 

By Lemma 8.2, it suffices to show (5.3) for the case that [L,L] = J^i^l^ %i 
is a compound-Poisson process, with jump times (r n ) ng jj. Then it follows 
easily by induction from (2.3) and (8.1) that 

N(t) 

Y t = e Bt Y +Y / e B ^eV ri Z l , t > 0. 

i=l 

In view of the proof of (b) below, let s > 0. Then 

N(t) 

(10.1) a!e Bs Y t = a'e^+OYo + £ a'e^+^Wr^ 

i=l 

N(t) 

(10.2) >7+E a ' eB(S+ ^ ri)e ^^. 

i=l 

Setting s = 0, it follows that Vt = ao + a'Y t _ > «o + 7 for t £ [0, Ti]; hence 
also Vr 1 + > «o + 7 > by (5.1) and (10.2), and an induction argument shows 
that Vt > ao + 7 for all t > 0, that is, (5.3) holds. 

For the converse, suppose first that (5.2) fails. Then, using the continuity 
of the function 1 1— ► e Bt , it follows that there is (£1,^2) C (0, 00) such that 
P(a + a'e st Y < Vt E (h,t 2 )) > 0, and since P(Ti > i 2 ) > we get the 
claim from (10.1). So suppose that (5.2) holds with 7 > — ao, but (5.1) fails. 
Suppose that the support of the Levy measure of the compound-Poisson 
process [L,L] (and hence the support distribution of the jumps, Zj) is un- 
bounded. Let (£3, £4) C (0,oo) be an interval such that a'e Bt e < — c\ < 
for all t £ (£3,24) for some c\ < 0. Let £5 > £4. By (5.2) we have P(Vr 1 > 
a o + 7) = 1) so that it is easy to see that the set 

A := {Ti < £5 < T 2 , t s - ri G (i 3 , i 4 ), Vr x > «o + 7} 

has positive probability. On A we have, by (10.1), 

V tB = a + a'e^Yo + a! e B ^~ Vl) gVt.Zl 
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Now a'e B ( t5 ~ ri ^e < —c\ and by choosing Z\ (which is independent of Ti,^ 
and Yo) large enough, we obtain P(Vt 5 < 0) > 0. 

(b) In view of (a) it remains to show that Yqo satisfies (5.2). For the proof 
of this, it suffices by Lemma 9.1 to assume that [L,L] is compound-Poisson. 
Let (Yj)t>o be a state process with Yo = 0. Then (5.2) holds for Yo with 
7 = 0, and it follows from (10.2), (5.1) and (5.3) that a'e Bs Y t > for all 
s,t > 0. Since Y t converges in distribution to Yoo as t — ► oo, (5.2) follows 
with 7 = 0. □ 



11. Proof for Section 6. 



Proof of Theorem 6.1. We mimic the proof of Proposition 5.1 of [19], 
that is, in the COGARCH(l,l) case. Observe that (6.1) and (6.3) follow 
immediately, since (Lt)t>o is a zero-mean martingale. Furthermore, (Gt)t>o 
is a square integrable martingale so that 



EG 2 r = E f V s d[L,L} s = E(L 1 ) 2 rE(V (x 
Jo 



and (6.2) follows from (4.9). Before showing (6.4), we verify that EG\ < oo 
if (4.3) is satisfied: it follows from the Burkholder-Davis-Gundy inequal- 
ity (see, e.g., [24], page 222) that EG\ < oo if E[G,G] 2 < oo. Let V t = 
a o|| e l|s,r + || ea/ ||_B,rYt_ be the volatility of the COGARCH(l,l) process 

constructed in Lemma 9.1 and let Gt = J * \fv~tdLt be the corresponding 
COGARCH(l, 1) price process. Then it follows from Lemma 9.1 that there 
is C\ > such that 

< V s = a + a'Y s _ < a + dY s _ = a + cj l 



l ea 'lb,- 



Then 



[G,G] t = fv s d[L,L\ 
Jo 



5 -nSh- l'v,^ L] , + U- ) [L, L]« 



ea'|| B ,rJ0 \ ll ea lis 



Cl -[G,G], + fa -^M^)[L,L] t , 



ll ea 'llB,r V ||ea'||s )T . 

so that again by the Burkholder-Davis-Gundy inequality and Doob's max- 
imal inequality, finiteness of EG t implies finiteness of E[G,G] 2 and hence 
of EGf. The fact that EG t < oo was used in [19] in the case when L has no 
Gaussian component, but it also holds in the general case. 
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Denote by E r the conditional expectation with respect to the a-algebra T T . 
Using partial integration, we have 

i-h+r 



{Gt ] ) 2 = 2f r G s .dG s + [G,G] h h X r 

Jh+ 

ph+r ph+r 

= 2 G s ^VV s dL s + V s d[L,L}t 

Jh Jh+ 



Since the increments of L on the interval (h, h + r] are independent of T r 
and since L has expectation 0, it follows that 

ph+r 

E r G s -W s dL s = 0. 

Jh+ 

Recall that Y s = J rjS Y r + K r s by (3.3) . Hence we also have Y s _ = J r s _ Y r + 
K T . jS _, so that, by the compensation formula, 

E r (G^f = E r [ k+r (a + a'Y,_)d[L,L] 

Jh+ 
rh+r 



rn+r 

= E r / (oo + a' J r , s _Y r + a'K rjS _) d[L, L] 

Jh+ 

(11.1) = E(Ll)a r + E{L\)el (EJ r>s _)Y r ds 

Jh+ 

+ E{L\)a! [ h+ \EK r ^)ds 

Jh+ 



ph+r 

E(Lj) / E r {V s )ds. 

Jh 



Since Yoo = J r ,s^x> + Kr,s by (3.5), with Yoo independent of ( J rjS , K r;S ) 

on the right-hand side, and EJ r ^ s = e B ( s ~ r ) by the proof of Theorem 4.2, it 
follows from (4.5) that 



Pq ~ «1M 



Hence 



E r (V s ) =a + z>e^Y r + e!-^—(I - e S ^) ei 

Pq ~ OL\\l 

(11.2) 

a °& +a ' e ^-)f Yr __^o^ ei 



/3 g - ai/i V Pq~ oi\ii 

Combining j£+ r e B(s-r) ds = e Bhg-i^_ e -B r ^ with (n.^ ( 1L2 ) and (4.5) 



gives 



E r {GW = E{L\) ( -p^- + a!e Bh B- 1 {I - e~ Br ){Y r - EY r ) 
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and we conclude with (6.2) that 
E{(Gt ] )\G^f) 

= E{E r {(Gl ) fGl)) 

= E(L\)e( a a ° rPq G 2 r + a 'e Sh B-\l - e- Sr )(Y r - EY r )G 2 r 
\p q - aifi 

= (E{G 2 .)) 2 + E(L 2 ) a 'e Sh B~ 1 (I- e~ Sr )[E(Y r G 2 ) - (EY r )(EG 2 r )], 
showing (6.4). □ 
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